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fT^ , Abstract. We have studied a system composed by two endohedral fullerene molecules. We have found that 

CO ' this system can be used as good candidate for the realization of Quantum Gates. Each of these molecules 

encapsules an atom carrying a spin, therefore they interact through the spin dipole interaction. We show 

!>■ 
O 



that a phase gate can be realized if we apply static and time dependent magnetic fields on each encased 
spin. We have evaluated the operational time of a vr-phase gate, which is of the order of ns. We made a 
comparison between the theoretical estimation of the gate time and the experimental decoherence time for 
each spin. The comparison shows that the spin relaxation time is much larger than the 7r-gate operational 



o 
o 



C ' time. Therefore, this indicates that, during the decoherence time, it is possible to perform some thousands 



^3 



of quantum computational operations. Moreover, through the study of concurrence, we get very good 
results for the entanglement degree of the two-qubit system. This finding opens a new avenue for the 
realization of Quantum Computers. 

PACS. PACS-key 03.67.-a ~ PACS-key 03.67.Lx - PACS-key 61.48.+C 

1 Introduction periments are aimed to realize a quantum gate. The first 

type of experiments is based on trapping ions in electro- 
During recent years there is a strong progress in model- ^.^gnetic traps, where the ions, which encode the qubit in 
ing physical realizations of a quantum computer. Many ^^^^^^ ^^^^^^^ freedom, are subjected to the mu- 
quantum physical systems have been investigated for the ^^^^ electrostatic interaction and to a state selective dis- 
realization of quantum gates. The most remarkable stud- placement generated by an external state dependent force 
ies were related to systems associated to Quantum Optics ^j^^ Cavity quantum electrodynamics (QED) tech- 
Ion Traps, to Quantum Electrodynamics in Optical Cav- ^.^^^^ ^^^^^ coherent interaction of a qubit, 
ities and to Nuclear Magnetic Resonance. All these ex- 
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generally represented by an atom or semiconductor dot 
system, with a single mode or a few modes of the elec- 
tromagnetic field inside a cavity. Depending on the par- 
ticular system, the qubit can be represented by the polar- 
ization states of a single photon or by two excited states 
of an atom. Although cavity QED experiments are very 
promising, they have been accomplished for few qubits 
[g]|^>|[7||§) . In the third experiment, nuclear spins repre- 
sent qubits. These spins can be manipulated using nuclear 
magnetic resonance techniques, and through the study of 
the quantum behavior of spins, quantum operations are 
realized. However, the number of spins which can be col- 
lected in a system is very limited, and this forbids the 
building up of a scalable quantum computer |^1(JII11II12I . 
From the study of such systems, we learn that the de- 
coherence phenomenon is the main issue which prevents 
the realization of quantum gates. Here we will focus on a 
physical systems, which will be able to produce a realis- 
tic quantum gate. The basic elements of our system are 
fuUerene molecules with encapsulated atoms or ions, which 
are called buckyballs or endohedral fuUerenes. Each of the 
trapped atoms carries a spin. This spin, associated with 
electronic degrees of freedom, encodes the qubit. It has 
been shown that these endohedral systems provide a 
long lifetime for the trapped spins and that the fuUerene 
molecules represent a good sheltering environment for the 
very sensible spins trapped inside. These endohedral sys- 
tems are typically characterized by two relaxation times. 
The first is Ti, which is due to the interactions between 
a spin and the surrounding environment. The second one 
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is T2 and it is due to the dipolar interaction between the 
qubit encoding spin and the surrounding endohedral spins 
randomly distributed in the sample. While Ti is dependent 
on temperature, T2 is practically independent of it. The 
experimental measure of the two relaxation times shows 
that Ti increases with decreasing temperature from about 
100/is at r = 3Q0K to several seconds below T = 5K, 
and that the value of the other relaxation time, T2, re- 
mains constant, that is T2 ~ 20/is |14II15| . In comparison 
with T2 the value of Ti is very large, therefore the system 
decoherence is determined by the spin-spin relaxation pro- 
cesses. It is supposed that the value of T2 can be increased, 
if it will be possible to design a careful experimental ar- 
chitecture, which could screen the interaction of the spins 
with the surrounding magnetic moments. It should be pos- 
sible to reduce the relaxation time of the system due to 
the random spin-spin interactions, if we consider a sys- 
tem composed by arrays of endohedrals encapsulated in 
a nanotube this system is also called as peapod, or 
considering buckyballs embedded on a substrate. These 
should be reliable systems for the realization of quantum 
gates. In such architectures the decoherence time for each 
encapsulated spin should be longer. 

Quantum computing through the study of doped fuUerene 
systems has been investigated in many works |17II18II19I 
1201121) . Although we have followed many ideas suggested 
in these previous papers, we consider a different approach 
for the realization of quantum gates. 

Our study is focused on a system composed by two 
buckyballs. Our aim is the realization of a quantum tt- 
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gate, which is a generahzation of the phase gate, this wiU 
be treated in Sec. To perform the 7r-gate, we need to 
know the time evolution of the coefficients of the standard 
computational basis states over which we expand the wave 
function of our system. The two particle phases are eval- 
uated through the numerical solution of the Schrodinger 
equation, see Sees. I5I6I We have used two approaches: a 
time independent Hamiltonian, see Sec. El and a time de- 
pendent one, see Sec. |B1 The main result of our study is 
the gate time, that is the time required by the system 
in order to perform the 7r-gate. The values obtained are 
around r ~ 1 x lO^^s, which is a few orders smaller than 
the shortest relaxation time, T2. From the comparison of 
the gate time, r, to the relaxation time, T2, we get that it 
is theoretically possible to realize some thousands of basic 
gate operations before the system decoheres. We have also 
checked the reliability of our gate through the analysis of 
the concurrence of the two-qubit state, see Sec. 0] The 
best value for the concurrence is obtained in the case of a 
time dependent Hamiltonian, while the gate time is nearly 
the same in both cases. 

2 Physical Features of the System 

The system under consideration is composed by two in- 
teracting buckyballs. Several experimental and theoreti- 
cal studies on buckyballs |13II17II22II23II24II25) . show that 
many different types of atoms can be encased in fuUerenes 
molecules. However, in most of the studied endohedral 
fuUerenes, there is a charge transfer from the encapsulated 
atom to the fuUerene cage, with a resulting considerable 
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alteration of the electronic properties of the cage. This 
is not the case for group V encased atoms. These atoms 
reside just at the center of the fuUerene molecule, there- 
fore there is no hybrididazion of the electron cloud of the 
encased atom and there is no Coulomb interaction with 
the fuUerene cage. In particular, the most promising en- 
dohedral molecule should be the N@Ceo, which is charac- 
terized by many interesting chemical-physical properties. 
Following Refs. |13II17II25| , experiments and theoretical 
calculations suggest that there is a repulsive exchange in- 
teraction between the fuUerene and the electronic cloud of 
the encapsulated atom. The electrons in the cloud of the 
encased nitrogen are tighter bound than in a free nitrogen 
atom, which allow the encased nitrogen to be less reac- 
tive even at room temperature. These results, together 
with the location of the nitrogen atom in the central site, 
suggest that in N@Cqo the nitrogen can be considered 
as an independent particle, with all the properties of the 
free atom. Since any charge interaction is screened, the 
fuUerene cage does not take any part in the interaction 
process and it can be considered just as a trap for the ni- 
trogen atom. Therefore, the only physical quantity of in- 
terest is the spin of the trapped particle. A nitrogen atom 
can be effectively described as a |-spin particle. This spin 
is associated with the electronic degrees of freedom. Tak- 
ing into account also the nuclear spin, which is ^ for the 
N@Cqq, the number of relevant degrees of freedom will be 
not increased [25] . We will consider a more simple model 
assuming that the encased atoms are described as ^-spin 
particles. In absence of any mutual interaction and with- 
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1 10) - 1 11); (4) 

I 11) -I 10). (5) 



out any applied magnetic field, the energy levels associ- 
ated with these spin particles are degenerate. If we apply 
a static magnetic field, this degeneracy is lifted. As a re- 
sult, due to the Zeeman effect, a two level system arises 
for each i-spin particle. Each of these two levels encodes 
the qubit. The spin- up component, nis = encodes the 
computational basis state | 1) , and the spin-down compo- 
nent, TTis — represents the state | 0). 

3 Gate Operation: The Phase Gate 

Quantum computers operate with the use of Quantum 
Gates. Quantum gates are defined as fundamental quan- 
tum computational operations. They are presented as uni- 
tary transformations, which act on the quantum states, 
which describe the qubits. Therefore a quantum computer 
must operate with the use of many quantum gates. The 
simplest gates are the single-qubit gates. Since our system 
is composed by two qubits, we will consider a two-qubit 
quantum gate. One of the most important quantum gates 
is the Universal Two-Qubit Quantum Gate 12 , which is 
called the CNOT-gate. The CNOT operation is defined 
by the following four by four unitary matrix 



'^1 0^ 



U. 



CNOT 



(1) 



10 
1 

yo 1 

and its action over the computational basis states reads: 



7 



00) 
01) 



00) ; 

01) ; 



(2) 
(3) 



The CNOT gate is given by the composition of a single- 
qubit Hadamard gate followed by a two-qubit 7r-gate, fi- 
nally followed by another single-qubit Hadamard gate. 
The representation of the Hadamard gate in the Bloch 
sphere is a ^ rotation about the y axis, followed by a re- 
flection of the X — y plane. In this paper we will focus 
on the realization of the two-qubit 7r-gate. It is a partic- 
ular choice of the general phase gate, represented by the 
following matrix 



G 



1 





1 





1 











(6) 



and its action on the computational basis states is the 
following: 



|00) ^ |00) 
101) - 101) 

|io) - |io) 



111) ^ e"'|ll) 



(7) 
(8) 
(9) 

(10) 



When d ~ ±7r, the resulting quantum gate is called a 
n~gate. In general, the time evolution of the four states 
of the standard computational basis can be described as 
follows: 



100)- 


-> e'-^^lOO) 


(11) 


|01)- 


e'"^"! |01) 


(12) 


|10)- 




(13) 
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(14) 4 Concurrence 



In order to obtain the action of the ideal quantum phase 
gate, equations (|7I10|I . see Ref. |3] , we have to apply the 
following local operator: 



(15) 



where 



^1 -I 0)1(0 I 



1)1(1 I 



S2 =\ 0)2(0 I e*^^+ I 1)2(1 I 
and the phases are defined as follows:: 

s? = ~(^oo/2 

Si = -'/'lo + 0oo/2 

s° = -4,00/2 

-(/)oi + 0oo/2, 



(16) 
(17) 



4 



(18) 
(19) 
(20) 
(21) 



After a straightforward calculation we obtain the desirable 
phase: 

■i^ = 011 - 010 - 001 + 000- (22) 

In our system, in order to realize a 7r-gate, we need to know 
the time evolution of the wave function. The time evolved 
wave function, expanded on the standard computational 
basis, is given by the following equation: 



iPit)) =ci(t) |00)+C2(i) |01) 

+C3{t) I 10) +C4(t) I 11). 



(23) 



Each coefficient Ci{t), i — 1, ..4, is a complex number, 
whose phase, arranged as in equation H22|l . is used for the 
realization of the 7r-gate. 



When we consider a ^-spin particle as the encoding system 
for the qubit, it may incur to a spin-flip process. This 
phenomenon consists in the swapping between the spin- 
up and spin-down components 



0) - I 1), 

1) -|0). 



(24) 
(25) 



If we consider the two-qubit state, known as EPR pair, 
I 00)+ I 11) 



V2 



(26) 



we can see that it is unaffected by the spin-flip of both 
qubits. This state, for this feature, is called maximally 
entangled. Therefore, we can define the entanglement as 
the property of quantum states, which shows if the state is 
good for carrying quantum information. The most entan- 
gled a quantum state is, the most reliable it is for transfer- 
ring quantum information. In our study we have consid- 
ered the concurrence, see Ref. |27) . as a measure of the en- 
tanglement of the state describing the two-qubit system. 
A pure state of two particles quantum system is called 
entangled if it cannot be factorisable, that is it cannot 
be written as the direct product of the states describing 
each particle. A mixed state if it cannot be represented as 
a mixture of factorisable pure states. In this Section we 
will refer to the entanglement of formation, which quanti- 
fies the resources needed for the creation of an entangled 
state. For a complete treatment about the entanglement 
of formation of pure and mixed states see Refs. |28II29| . 
The entanglement of formation of a quantum state can 
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be evaluated through the concurrence |27| . Since the state Finally, we obtain the concurrence of our system, see eq. 

describing our system is a pure state, the degree of en- (|27|l . by performing the state product between states (|32|l 

tanglement of our system can be quantified through the and The normalized concurrence of the system is 

definition of the concurrence for a pure state j27| . which given by the following equation 

is defined by o\ r*r* r*r* I 

Ci^) = , „ , ^ f ' , „ . (34) 

C(^) =1 I ^) I, (27) |ciP + |c3P + |c3P + |c4P 

, i7\-,i -a- jixr i rni -a- The rcsult obtained in eq. (I34II will be used to evaluate 
where | %!) is the spm-nipped state oi system. Ihe spm-nip 

r i-ir 1 • ijj the degree of entanglement of our system during the gate 

transformation, which for a ^-spm particle is the standard ° ° j a o 

, , „ , . • 1 r. 1 r 11 operation. When the concurrence related to a wave func- 

time reversal transformation is denned as follows 

tion reaches its maximum value, the state is maximally 

\i')^ay\ r), (28) 

entangled. Therefore, we require that the concurrence of 
where CTy is the Pauli y-matrix and | ip*) is the complex the wave function of the system, at the end of the gate 
conjugate of | tA)- The entanglement, see EZj, is defined as operation, reaches a value next to its maximum, 
a function of concurrence, through the following equation 

E{^j) = f{C{^)), (29) 5 Phase Gate: Time Independent Case 

where function fiCiib)) is given by ^ ^ 

^ ^^^^ * 5.1 Preliminary Setup 

Our system is composed by two spins, which interact with 



= + ,30) 



h(x) = —x\og2X — (1 — x)\og2{l — x). (31) 

a static magnetic field. Applying a static magnetic field 

Function f{C{ip)) increases monotonically from to 1 as oriented in the z direction, for the Zeeman effect, we get 

C{ip) ranges from to 1. Therefore, the concurrence can the splitting of the spin z component into the spin-up and 

be considered as a measure of the entanglement. spin-down components. The energy difference between the 

The state describing our two-qubit system, written as a su- two levels give the resonance frequency of the particle. 

perposition of the standard two-qubit computational basis However, when we apply a static magnetic field on the 

states, is given by whole sample, all the particles will have the same reso- 

, , , , , , , , , , , , ^ nance frequency. To perform manipulations on each buck- 
|V')=ci|00)+C2|01)+C3|10)+C4|ll). (32) H y y 

yball, we need to be able to distinguish each of them. 

Following eq. 1)2811 . the spin-flip transformation over the 

This setup leads to the most relevant experimental dis- 

state gives 

advantage for systems composed by arrays of buckyballs, 
I ■0) = —cl I 00) -I- C2 I 01) + Cg I 10) — C4 I 11). (33) which is the difhculty in the individual addressing of each 
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qubit particle. This problem can be overcome with the use 
of external field gradients, which can shift the electronic 
resonance frequency of the qubit-encoding spins |18II2U) . 
Magnetic field gradients can be generated by considering 
wires through which flows current. If we place two parallel 
wires outside our two buckyball system, it is generated an 
additional magnetic field in the space between the wires. 
Following a paper by Groth et al. (211, with the help of 
atom chip technology, wires with a high current density 
can be built. The magnetic field amplitude generated by 
the two wires is given by 

^'^2i^^^x + p + d/2^ x-p- ^^^^ 

where / is the current intensity, d is the distance between 
the two wires, p is the radius of each wire and x is the 
distance of a buckyball with respect to the origin of the 
axes. With the choice / = 0.6yl, d — 1pm and p = l/im, 
through a numerical computation, we obtain the magnetic 
field distribution shown in Fig. We could not consider 
a current greater than / = 0.6A because the wires would 
face a too high heating process, and eventually they could 
be destroyed. On the other hand, we could not consider 
currents smaller than lO^^A, because the arising mag- 
netic field gradient would be too small for each buckyball. 
In this case, the resonance frequencies related to the buck- 
yballs would differ for only few Ad Hz, which could be a 
too small gap to be realized by a frequency resonator. 
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Fig. 1. Magnetic field generated by two 1/im-radius wires at 
a distance d = 1pm, which carry a current I — 0.6A. The two 
buckyballs are placed at a symmetrical distance x with respect 
to the origin of the axes. 

5.2 Realization of the Phase Gate 

Choosing a static magnetic field in the z direction, the 
Hamiltonian of the system is given by the following equa- 
tion {h — 1) 

H = Jo&i ■ (T2 + g(r)[(Ti • (T2 — 3((Ti • n){a-2 ■ n)] 

-PB[{{B,,+Bg,)(7,,)®h (36) 

In the previous equation, Jq is the exchange spin-spin in- 
teraction coupling constant, cti and (72 are the Pauli spin 

2 

matrices, g{r) — 7172^7^, where /io is the diamagnetic 
constant, pB is the Bohr magneton and r is the distance 
between the two trapped atoms, n is the unit vector in 
the direction of the line which joins the centers of the two 
encased atoms, B^^ = B^^ is the static magnetic field in 
the z direction, Bg-^ and Bg^ are the additional magnetic 
fields due to the field gradient. We make an assumption, 
considering the trapped particles as electrons. Therefore 
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the gyromagnetic ratio 7 ~ 2, and g{r) — ^"^^f . Through 
the study of fuUerenes' spectra in ESR (Electron Spin Res- 
onance) experiments, and also through theoretical studies, 
it has been shown |13II32II17| . that the exchange interac- 
tion is very small. Therefore, in eq. (|36|l . we can neglect 
the exchange term proportional to Jq, leaving the spin 
dipole-dipole interaction as the leading term of the mu- 
tual interaction between the two endohedrals. Choosing 
the direction of vector n parallel to the x axis, the dipole- 
dipole interaction term is simplified as follows 



C2 W = -i[(-.9(0 + m2)c2{t) - .g(r)c3(i)]; (43) 
C3{t) = -i[-g{r)c2{t) + (-g(r) - m2)c3{t)]- (44) 
C4{t) = -i[-3gir)ci{t) + [g{r) - rm)ci{t)], (45) 

which allows us to evaluate the phases acquired by each 
computational basis state during the time evolution. Ap- 
plying eq. H22() to the present time evolved phases, we get 
the desirable 7r-gate 



matrix 



d = Arg(c^{t)) - Arg{c2(t)) (46) 
- Arg{c3{t)) + Arg{c^{t)) = ±7r, 

^ = 9{'^)(.^zi(^z2 + ^yi^yi ^ 2(5'a;i(7j;2). (37) where Arg{ci{t)), i — 1, ..,4, which correspond to phases 

4>ju j,l — 0, 1, in eq. H22|) . are the phases of coefficients 
Ci{t) of equation H23|) . We have numerically solved the dif- 
ferential equation system (|42I45|) , with the use of a Mathe- 
matica programme. The numerical quantities used for the 
numerical calculations are r = 1.14nm, Bzi ~ Bz2 — 



The Hamiltonian matrix form is given by the following 



g{r) + mi 

—gif) + 1^2 





-gir) 



-3gir) 




~g{r) ~gir) -m2 

-3g{r) g{r) - mi 



where 

mi = -usiBzi + Bg^ + Bz2 + SgJ 

and 

m2 = -fJ.B{Bzi + Bg^ - Bz2 - Bg.^), 

are the static magnetic field terms. Solving the 
Schrodinger equation 



where the wave function is a superposition of the standard 
two-qubit computational basis, given by equation (|23|l . we 
get the four differential equation system 

ci(i) = -iMr) + mi)ci{t) - 3.g(r)c4(t)]; 



lOx 10"2y, Bg^ = 6.08X lO^^T and Bg^ = -6.08x lO^^T, 

(38) 

which give the resonance frequencies tJi — 1.7599x 

and L02 — 1.7577X IQ^^Hz. The time evolution of the phase 

(39) 

•d is shown in Fig. |2] The gate time, which corresponds to 
the case i? = — tt is r ~ 9.1 x 10~®s. This result has been 
(40) found for a chosen set of initial conditions Ci(0), i = 1, .., 4. 
However, we did many trials for different numerical val- 
ues of the set Ci(0), i = 1, ..,4. In all these cases, phase 9 
shows a linear behavior and the resulting gate times are 
all in the same range, which is of the order of lO^^s. If the 
set of initial conditions is real, the starting value of phase 
9 is always equal to zero. If the set of initial conditions is 
complex, the starting value of 9 is in the range [— tt, -l-7r], 
(42) but it can always be rescaled to zero. The numerical value 



(41) 
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time (ns) 

Fig. 2. Time evolution of the phase i9. The value j9 = —vr is 
reached at the time t ~ 9.1 x 10~^s. 

of the distance between the two buckyballs, r, is a fixed 
value, which depends on the substrate where the buck- 
yballs reside. The amplitude of the static magnetic field 
has been found by considering the allowed experimental 
limits for its realization. The chosen value for this am- 
plitude has been found by checking the response of the 
system, i.e. the gate time, after some trials. Therefore, we 
can say that the phase gate time depends on the distance 
between the two buckyballs and on the amplitude of the 
static magnetic field, but it is independent of the choice 
of the initial values Ci(0), i = 1, ..4. 

If we compare the gate-time, r, to the shortest decoher- 
ence time, T2 — 20/iTO, we can deduce that it will be 
theoretically possible to realize about thousands gate op- 
erations before the system relaxes. To know the fidelity 
of the gate and the reliability of the results, we need to 
evaluate the concurrence during the time evolution. With 
the use of a Mathematica programme we have plotted 
the time evolution of the concurrence, equation H34|) . from 
t = Os to the gate time t = t, see Fig. O Analyzing pic- 
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Fig. 3. Time evolution of the concurrence, C{ip). 

ture |j3J|, we can see that the concurrence shows a smooth 
behavior. It monotonically ranges from zero and its max- 
imum is reached at time t = t, with the respective value 
C(V'(t)) = 0.88. Even if the maximum concurrence does 
not coincides with the ideal value 1 , it is near to this value 
and the system shows an acceptable degree of entangle- 
ment. It is convenient to investigate other system config- 
urations, in order to check if it is possible to improve the 
concurrence. In the next Section we will analyze the case 
of an additional magnetic field, oscillating in time in the 
x-y plane. 

6 Phase gate: Time Dependent Case. 

In this Section, we apply to our system an additional time 
dependent magnetic field. To induce the transitions be- 
tween the two Zeeman energy levels, we need to apply an 
oscillating magnetic field in the x ~ y plane with angular 
frequency, uj, equal to the spin resonance frequency. In the 
case of a transverse linear oscillating magnetic field, the 
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0. 3. 6. 9. 




0. 3. 6. 9. 



time {ns) 

Fig. 4. Time evolution of the phase with the respective 
gate time r ~ 9.8 x W~^s. 

total applied magnetic field is given by 

B(t) = {Bicosujt,Bicosujt,{B^ + Bg)). (47) 

The Hamiltonian of the system reads 

H = g{r){az^a2^ + Uy^ay^ - 2(7^^(7^^) 
-f^B{Bzi + Bgja:,, «) h 

-^iB{B,,+Bg,)h®a,, (48) 
-^bBi^ (cr^i COS Wit + Uy^ coswii) ® h 

+h ® {~^^BBl^{aj:^ C0SUJ2t + Oy.^ C0SW2t))- 

Like in the time independent case, solving the Schrodinger 
equation, we get a four differential equation system, whose 
solution give the time evolution of the phase for each com- 
putational basis state. Arranging the phases as prescribed 
in equation H22|) . we have obtained the 7r-gate. In the nu- 
merical computation we have used the additional quantity 
Bi^ = Bi^ = 5 X lO^^'r. shown in Fig. 21 and the numeri- 
cal value of the gate time is r ~ 9.8 x lO^^s. Also in this 
case, comparing the gate time, r, to the decoherence time 
T2, we observe that it will be possible to perform about 
thousands gate operations before the system relaxes. The 
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Fig. 5. Time evolution of the concurrence, C{ip). 

relevant result in the treatment of the time dependent case 
is the concurrence. In Fig. |31 it is represented the time 
evolution of the concurrence, C('(/'(t)), which has been nu- 
merically evaluated with a Mathematica programme. It 
shows a monotonic behavior and its maximum, evaluated 
at time t — t, corresponds to C('0(r)) = 0.96. There- 
fore, an additional linearly polarized oscillating field in 
the X — y plane allows the system to be characterized by 
a better concurrence degree. 

7 Conclusions 

To model quantum gates we considered a system com- 
posed by two endohedral fuUerene molecules, subjected to 
external magnetic fields. We assume that each molecule 
may be treated as a i-spin particle, where the spin is as- 
sociated to the encapsulated atoms. In the magnetic field 
the spin degeneracy of the spin up and down components 
is lifted and it arises the Zeeman splitting. As the result, 
there two two-level system are arising. Each of these two- 
level systems corresponds to a single qubit. If the applied 
static magnetic field to the whole sample is homogeneous. 
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each of these qubits will be characterized by the same 
resonance frequency. This leads to the difficulty in the in- 
dividual addressing of each single qubit. To overcome this 
problem, we have to apply inhomogeneous magnetic fields, 
in this paper we have used a magnetic field generated by 
two metallic wires. Each wire is carrying a current, there- 
fore the magnetic field is decreasing with the distance from 
a wire. In the proposed configuration of two parallel wires, 
there arises a gradient of the magnetic field when we are 
moving from a wire to the other one. If we place two buck- 
yballs in the space between these two wires, they will be 
subjected to the gradient of this field, and therefore the 
associated resonance frequencies of the related two-level 
system are different. In this paper we have performed a 
quantum vr-phase gate. To realize this particular quantum 
gate we have estimated the phase of each computational 
basis state, see equation (|22|) . The leading mutual inter- 
action between the two qubits is the spin dipole-dipole 
interaction. First we studied the time evolution of our sys- 
tem taking into account this mutual interaction between 
the qubits and considering the qubits subjected to static 
magnetic fields only. Then we applied to the system also 
time dependent magnetic fields. The wave function of the 
system is given by the superposition of the four computa- 
tional basis states, see equation (|23|) . The time evolution of 
the coefficients of each computational basis state is deter- 
mined via the solution of the Schrodinger equation. With 
the use of these coefficients and of equation (|46|) , we can 
evaluate the operational gate time for the 7r-phase gate. 
Its numerical value is r ~ 9.1 x lO^^s for the time inde- 
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pendent case, and t ~ 9.8 x lO^^s for the time dependent 
one. Comparing the gate time, r, to the shortest relaxation 
time, T2, we have observed that in both cases it will be 
possible to perform about thousands quantum gate opera- 
tions before the system decoheres. This is our main result. 
As far as we are aware, this result indicates that our sys- 
tem could be the most favorable for the realization of a 
quantum gate. Obviously, for realistic models of quantum 
computers, the ratio of the decoherence time and the op- 
erational time must be very large, otherwise the system 
relaxes before the completing of the quantum computa- 
tion. The goal of any quantum computational proposal is 
the entanglement of the state of the system under consid- 
eration. At this purpose, we have studied the concurrence, 
see Sec. 0] The concurrence gives information about the 
entanglement of the state, therefore it is related to the 
reliability of the gate operation. A maximally entangled 
state is left unchanged under a spin-flip operation and its 
concurrence is maximum. In our system, at the end of the 
gate operation, the value of the concurrence is C ~ 0.88 
in the time independent case, and C ~ 0.96 in the time 
dependent one. Both values are acceptable because they 
are both related to a very good degree of entanglement 
for the state describing our system. We can conclude that 
the best configuration for our system is the time depen- 
dent one. It is characterized by a very small operational 
time, in comparison to the relaxation times, and by the 
best concurrence. 

Many features claim the buckyball systems as good can- 
didates for performing quantum gates. Not only they are 



12 Maria Silvia Garelli, Feodor V Kusmartsev: Buckyball 

characterized by very long decoherence times, but also 
they can be maneuvered very easily. This feature allows 
the realization of experimental quantum devices, which 
form scalable architectures. For example, buckyballs can 
be embedded in silicon surfaces or arranged in arrays en- 
cased in a nanotube (peapod). Moreover, in such systems 
we suppose tliat the value of the relaxation time T2, duo to 
random spin dipole-dipole interactions, could be reduced. 
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